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Abstract. The properties of materials near structural phase transitions are 
often successfully described in the framework of Landau theory. While usually 
the focus is on phase transitions, which are induced by temperature changes 
approaching a critical temperature Tc, here we will discuss structural phase 
transformations driven by high hydrostatic pressure, as they are of major 
importance for understanding processes in the interior of our earth. Since at 
very high pressures the deformations of a material are generally very large, 
one needs to apply a fully nonlinear description taking physical as well as 
geometrical nonlinearities (finite strains) into account. In particular it is necessary 
to retune conventional Landau theory to describe such phase transitions. In 
[A. Troster, W. Schranz and R. Miletich, Phys. Rev. Lett.88, 55503 (2002)] we 
have constructed a Landau-type free energy based on an order parameter part, an 
order parameter- (finite)strain coupling and a nonlinear elastic term. This model 
provides an excellent and efficient framework for the systematic study of phase 
transformations for a wide range of materials up to ultrahigh pressures. 
We illustrate the model on the specific example of BaCr(Si40io), showing that it 
fully accounts for the elastic softening which is observed near the pressure induced 
phase transformation. 
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1. Introduction 

A major achievement of Ekhard Salje and coworkers was to show that Landau theory 
can be very efficiently applied for the desription of experimental data near temperature 
induced structural phase transitions in minerals (Salje, 1990; Salje, 1992; Carpenter, 
Salje and Graeme-Barber, 1998). There are a couple of reasons for this success. First 
of all, as coupling to strains induce long range forces and enhance anisotropy, strain- 
induced interactions always tend to suppress fluctuations, leading to a mean field type 
of transition (Bratkovsky, 1995, Tsatskis, 1994). Depending on the nature of the 
system (i.e. on the type of coupling between strain and order parameter its anisotropy 
and the boundary conditions imposed on its surface) strain effects also may change 
the character of the transition from second to first order (Bergman and Halperin, 
1976). Now in principle, a Landau potential can be obtained from a low order Taylor 
expansion of the mean field free energy. The problem of how to determine the range 
of validity of this expansion naturally arises. To answer this question one observes 
that in many mineral systems the phase transitions turn out to be of the displacive 
type (Dove, 1997). For such systems, the above authors were able to show that 
Landau theory indeed allows to reproduce experimental data in an extremely broad 
temperature range, in particular down to very low temperatures by taking into account 
the effect of quantum saturation (Salje, Wruck and Thomas, 1991). The situation is 
not so simple for the other extreme, i.e. order-disorder phase transitions, since it is 
no longer possible to approximate their free energy by a low order polynomial in the 
order parameter (Giddy, Dove and Heine, 1989). In passing, we note that real systems 
are often of mixed, i.e. displacive/order-disorder character (Meyer, 2000; Sondcrgcld, 
2000), and the characterization of crossover systems in between a displacive and order- 
disorder type is still a matter of active research (Perez-Mato, 2000; Rubtsov, Hlinka 
and Janssen, 2000; Troster, Dellago and Schranz, 2005). 

In "traditional" Landau theory it is difficult to include pressure effects beyond the 
infinitesimal strain approximation. Therefore one frequently resorts to replacing the 
elastic background energy by its harmonic truncation and consequently treats the 
strains as infinitesimal. Since the total strains appearing in temperature induced 
phase transitions arc as a rule quite small this is well justified. However, when we 
discuss the role of extremely high pressure in driving the transition, the validity of 
this assumption must be seriously doubted. 

At high pressure a phase transition is usually detected by anomalies e. g. in the 
system volume V(P), which in turn result from anomalies in the pressure dependence 
of the lattice parameters ai{P), i = 1,2, 3, near a critical pressure Pc- In many cases 
this leads to an additional (with respect to the " background" ) nonlinear P-dependence 
of lattice parameters and/or volume. As a function of pressure, the "stiffness" of 
any solid is characterized by the isothermal compressibility k{P) = —d\ogV{P)/dP. 
Various theoretical concepts are usually employed to derive so-called equations of state 
(EOS) (Anderson, 1995; Angel, 2000a; Angel, 2000b), which - in absence of phase 
transitions - describe the hydrostatic pressure dependence of the crystal's volume 
V{P). 

In presence of a high pressure phase transition (HPPT) a commonly adopted 
practice is to merely fit the corresponding V'(P)-behavior to a number of differently 
parametrized EsOS for each phase (Schulte and Holzapfel, 1995; Kriiger and Holzapfel, 
1992; Chesnut and Vohra, 2000), and, apart from computer simulations, in the case of 
a reconstructive phase transition there is little more one can do at present. However, 
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from a theoretical point of view, although this purely phenomenological procedure 
does fit experimental volume data in many cases, it neither describes the pressure 
behaviour of individual strain components nor provides any possible further insight 
into the mechanisms of a HPPT. A more profound theoretical approach is therefore of 
vital interest to a broad audience reaching form physicists studying the high pressure 
behavior of materials (crystals, liquid crystals, complex liquids, biological membranes, 
etc.) to geologists investigating the earth's minerals and bulk properties. 

Indeed, for HPPT's of the group-subgroup type, it is natural to expect that one 
could do much better by applying the concepts of Landau theory. In particular, 
we should be able to connect the high and low pressure phase by a uniform 
thermodynamic description. However, here one faces a marked difference compared 
to the temperature driven case. At high pressures, as the inter-atomic forces opposing 
further compression, a crystal's volume and lattice parameters develop large strains 
in a pronounced non-linear way, which corrupts any serious attempt to treat the 
elastic energy in the framework of the infinitesimal strain approximation. For these 
reasons it is necessary to reconsider the construction of Landau theory at extremely 
high pressures and large deformations (Troster, Schranz and Miletich, 2002). Here 
we summarize the main ideas of this new theoretical approach and discuss some 
consequences for the elastic behaviour of materials. 

2. New approach for high pressure phase transformations 

As stated above, the conventional infinitesimal strain approach is clearly insufficient 
for aiming at an accurate description of HPPT's. Instead, at high pressures the elastic 
response of a system must be characterized in terms of a appropriate nonlinear strain 
measure like the Lagrangian or Eulerian strain tensor. Let X and X' denote an 
undeformed reference system and a deformed system, respectively. Then, from the 
deformation gradient tensor aik = dX[/dXk one constructs the Lagrangian strain 
tensor eik = \ (Sn '^'"'-'^nfc ~ ^ik)- In the present work, we will consider three pairs 
of systems X, X' at hydrostatic external pressures P (FigH]). For measuring the 
"background strain" e(P), which is defined as the Lagrangian strain displayed by the 
system at the constraint of zero order parameter Q (see below), X{P) denotes the 
corresponding state and X{P = 0) is used as a reference system. For measuring 
the spontaneous strain e resulting from the emergence of Q 7^ 0, X{P) plays the 
role of a (hypothetical) "floating" background reference system. On the other hand, 
above the critical pressure Pc the total strain r]{P), which is obtained as the nonlinear 
superposition rj = e + a~^.e.a. of the spontaneous strain e appearing "on top of" 
the background strain e(P), is the strain measured in the fully deformed state Y{P), 
which is the one actually displayed by the system, the system X{0) again playing the 
role of the reference system. 

In the spirit of Landau's theory of phase transitions the Landau free energy 
F{Q,e) is built from an order parameter part Fq{Q), an order parameter-strain 
coupling FQe{Q, e; X{P)) and a purely elastic free energy term Fe(e; X{P)) (Toledano 
and Toledano, 1988). The order parameter part Fq{Q) is usually written as a 
polynomial in the variables Q which is constructed using the theory of invariants of 
irreducible representations of the space groups (Kovalev, 1993). The construction of 
the coupling term Fq^iQ, e; X{P)) is more tricky. Measuring the strains with respect 
to the floating background reference system X{P) and assuming linear-quadratic 
coupling of strain and order parameter with coupling coefficients dlHX{P)), we write: 
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Figure 1. Sketch of deformation states and superposition of strains. 



FiQ, e; X(P)) = T/(P)<I'(Q; X(P)) + ^ 4^(X(P))Q,Qj% 
+ Fo(e;X(P)) (1) 



HPPT's are characterized by the fact that the total strain r] is not small. However, 
for phase transitions close to second order there is at least a nonzero pressure interval 
starting at Pc where the spontaneous strain e induced by the nonzero equilibrium 
order parameter Q can still be treated as infinitesimal. In this approximation the 
pure elastic free energy part of the Landau potential with respect to X{P) reads: 



Fo{e; X(P)) « T/(P) j ^ r,,%i ^ aM^iPm, 

\ ij ijkl 



(2) 



where Cijki{X{P)) are the crystal's thermodynamic elastic constants at pressure 
P. Inserting Eql2]into Eqlljand minimizing with respect to eij yields: 

Y,df^HX{P))QKQL+Y.aMX{P))eki{Q) = (3) 

KL kl 

Solving Eql3] for Cij (P) and inserting into the equilibrium equation for the order 
parameter 
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- ^''^y^'^^^ + 2 5:4-(X(P))Q,?,(Q) (4) 

one obtains the renormalized order parameter potential written with respect to 
the deformed state X.[PY 



$fl,(Q;X(P)) :=$(Q;X(P))- 

-\Y. QiQjQkQl [Y^4';{X{P))C-^l{P)d^i\X{P))\ (5) 

IJKL \ijkl J 

In the spirit of Landau theory we now make the simple but crucial assumption 
that the Landau and coupling coefficients of the potential depend on P only in a 
"geometrical" way, i.e. through volume ratios and geometrical transformation rules 
between the reference states X{P) and X{0), such that the P dependence of the 
equlibrium order pararameter is due to the presence of the order parameter-strain 
coupling. Transforming Eq|S]to the undeformed reference state X(0), the equilibrium 
order parameter Q can then be calculated as the minimum of 

<^>r{Q;X{0)) = ^Q;X{Oj) + Y,dlf {X{Q))e,,{P)QiQj - 

ijij 

- dif{Xm'^^^^^d^i^{X{0))QjQjQKQL (6) 



IJKL 

ijkl 



where 



Trost{X{P)) = yfm X! '^^rajoC^jli{X{P))aksait (7a) 

diiiXiP)) ^^^a™d,(/(X(0))a„, (76) 

ij 

and the total (Lagrange) strain r]ij turns out to be 



V^J = -Y.QKQLY.d^i'^iXiomMP) (8) 

KL kl 

Applying the above procedure and employing the usual assumptions concerning 
the temperature behaviour of the Landau parameters, one is able to calculate 
pressure and temperature dependencies of thermodynamic quantities like specific 
heat, strains and elastic constants, soft mode behaviour, etc. near HPPTs. A 
concrete example of the above theory at work will be given below. However, the 
above constructions obviously rely on the knowledge of the pressure dependence of 
the elastic constants Cijki{X{P)) = Cijki{P) (or equivalently, the Birch coefficients 
Bijki (P) = P{5ij5ki - SiiSjk - SikSji) + Cijki (P)) and the resulting background strains 
Bij (P) . Note that these quantities appear also in the renormalized free energy density 
Eql6l To calculate it we proceed in the following way. 
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For crystal classes that develop no shear strains under hydrostatic pressure (cubic, 
tetragonal, hexagonal, orthorhombic) the deformation tensor aij{P) = ai{P)6ij is 
diagonal and the axial compressibilities satisfy the equations 



with the initial conditions Q;i(0) — 1, where Sij{P) := B~^{P) denote the 
compliance tensor elements and we switched to Voigt notation. Therefore the 
deformation tensor components ai{P) as well as the finite strain tensor components 
eij{P) or lattice parameters ai{P) can be easily calulated by integrating EqlHl once 
the functions S'ifc(P), i,k = 1, 2, 3 are known. 

In our recent work (Troster, Schranz and Miletich, 2002) we proposed the 
expansion 



where k{P) is the bulk compressibility, kq = k{P = 0) and 5° denotes the zero- 
pressure compliance. The expansion coefficients obey the sum rule ^j^- = 
V n e N. Eq llOi which is exact at infinite n, yields an elegant parametrization 
of the pressure dependence of the compliance tensor, since it essentially factorizes 
out the main nonlinearity of the compliance through the bulk compressibility k{P)- 
The elastic anisotropy is taken into account by expanding the remainder in powers of 
P, as we anticipate that the very nature of the expansion allows it to be truncated 
at low orders with good accuracy. Indeed we have shown in a very recent work 
(Koppensteiner, Troster and Schranz, 2006) that even a low order truncation of our 
ansatz Ea llOl excellently fits the experimental high pressure data of olivine, fiuorite, 
garnet, magnesium oxide and stishovite. In fact, for all these examples a truncation 
at order n — 1 turned out to be sufficient to produce accurate results, the only 
exception being magnesium oxide, where the expansion Ea llOl had to be taken to 
n = 2. This very encouraging result implies that our new thermodynamic theory 
indeed allows for a dramatic reduction of the number of fit parameters as compared to 
conventional nonlinear approaches. Tablc[T] shows the numbers of relevant (nonshear) 
elastic constants of second, third and fourth order for various point group symmetries. 
Our approach requires only (72 — 1 additional unknown parameters njj (at order n = 1), 
whereas a fit of an expansion of the elastic energy including nonlinear elastic constants 
up to fourth order introduces 93 + additional fit parameters. For example, in the 
case of cubic I symmetry, our new approach requires 52 ^ 1 = 1 ht parameters as 
compared to 52 + 93 = 7 unknown parameters of a conventional fourth order nonlinear 
theory. For lower symmetry the advantages of our approach are even more obvious. 
For instance, in the case of orthorhombic symmetry our new parametrization needs 
5 parameters, which should be compared to a total of 25 unknowns for conventional 
nonlinear elasticity! 

The combination of high pressure Landau theory and the use of the expansion 
(jlOp for the background system X{P) is illustrated by an analysis of measurements 
of tetragonal BaCr(Si40io) single crystals. Using x-ray diffraction, the P-dependence 
of lattice parameters ai{P) — a2{P), a3{P), and the unit cell volumes V{P) were 
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Table 1. Numbers 92-93,54 of nonshear elastic constants of second, third and 
fourth order (Roy and Dasgupta, 1988; Brendel, 1979). 
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Figure 2. Pressure dependence of lattice parameters and unit cell volume of 
BaCr(Si40io). Points are measured data, lines are fits using the present theory. 

measured very detailed at room temperature in a diamond anvil cell (Troster, Schranz 
and Miletich, 2002). 

One finds a tetragonal-tetragonal HPPT at approximately 2.24 GPa, character- 
ized by a discontinuitiy in ai{P), 03 (P) and V{P). Consistent with the observed 
pressure hysteresis behavior, the transition can be classified as being of (weakly) first 
order. The order parameter part is constructed in the following standard way (Ko- 
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Table 2. Fit parameters for the high pressure Landau theory of BaCr{Si40io)- 



BaCr(SiiOio) 



A 0.45 GPa 

B -0.2 GPa 

C 20 GPa 
ai(0) = a2(0) 7.535 A 

03(0) 16.09 A 

K? 0.0035 GPa-i 

k9 0.001 GPa-i 

S\ 0.0035 GPa-i 

S§ 0.0089 GPa-i 

K 4.1 

41 1.9x10-"' GPa-2 

42 -1.7x10-'' GPa-2 
2.3x10-4 GPa-2 

di 13.75 GPa-i 

d3 -0.33 GPa-i 



valev, 1993): The symmetry reduction P4/ncc to P42i2 is driven by the onedimen- 
sional irreducible representation T2 at the wavevector fc = 0, yielding a one component 
order parameter Q, which is zero in the paraphase {P < Pc) and nonzero in the dis- 
torted phase (P > Pc). For ^{Q,X) we assume 

where A, C > 0. The tetragonal symmetry also dictates di{X) — d2{X) ^ 

Let Kq := K~^{0) denote the isothermal bulk modulus at P = 0. The Murnaghan 
equation of state (MEOS) (Anderson, 1995) 

viP) ■■-y^ = (^ + KP/Ko)-'"''" (12) 

which is based on the simple ansatz k^^{P) —: K{P) = Kq + KqP, is frequently 
used to describe (P, y)-data and is known to usually reproduce the values of K{P) 
correctly up to volume changes somewhat larger than w(P) > 0.9 while being 
algebraically much simpler than other approaches such as the "Vinet" or the "Birch- 
Murnaghan" EsOS (Anderson, 1995; Angel, 2000b) used for higher compression 
ranges. Figl2] shows corresponding fits of unit cell volume and axes of BaCr(Si40io) 
using the parameter values of TablelSl Note, however, that our theory can in principle 
be used in combination with any of these EsOS, or even a function V{P) derived from 
an experimental measurement or a computer simulation. With these values our model 
confines possible pressure ranges for hysteresis effects to 2.2 GPa-2.4 GPa. One also 
calculates that the geometrical error introduced in assuming the spontaneous strain e 
to be infinitesimal is smaller than 0.9%, yielding an error < 0.1% in the total strain 
V- 

Figure [3] shows the pressure dependence of the bulk modulus of BaCr{SiiOiQ) 
calculated from the experimental data of V{P) as compared to the predictions of 
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Figure 3. Pressure dependence of bulk modulus K{P) calculated from the 
experimental data of V{P) (points). The lines show calculations based on 
piecewise EOS fitting (thin line) and the present high pressure Landau theory 
(thick line). 



the different approaches discussed above. It is clearly evident that the conventional 
fitting procedure, which is based on a piecewise fitting of EOS in the high and low 
pressure phase, respectively , underestimates the elastic anomaly by at least a factor of 
two, whereas the present high pressure adapted Landau theory reproduces the elastic 
anomaly very well. The same behaviour is found for the axial imcompressibilities 
(FigHl). It is worth noting that a very similar softening was also observed in other 
examples of pressure-induced phase transitions, notably for the bulk modulus of solid 
Ceo near its fcc-sc transition (Pintschovius, 1999) and for the longitudinal acoustic 
modes near the cubic-tetragonal transition of BaTiO^ (Ishidate, 1989). 

3. Summary 

Summarizing, our theory allows to compute a wealth of experimental observables 
(e.g. lattice parameters, elastic constants, specific heat, soft modes, etc.) from a 
quite transparent thermodynamic model based on coupling finite strain elasticity to 
Landau theory. There are several advantages of this new approach to HPPTs: At first 
there is a drastic reduction of fit parameters as compared to an expansion in terms 
of nonlinear elastic constants (Tabl^. Moreover the present approach allows for a 
direct connection to EOS fitting procedures. Compared to the method of piecwise 
EOS fitting the main advantage of our method is that it allows for a unified consistent 
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BaCrSip^ 




Figure 4. Pressure dependence of axial incompressibilites Ka , Kc calculated 
from the lattice parameters ai , 03 (points) and the bulk modulus K compared 
with the results of the present Landau theory (lines). 



description of high and low symmetry phases with a single set of parameters. 
At present, our theory is capable of dealing with group-subgroup transitions involving 
symmetries of cubic, tetragonal, orthorhombic or hexagonal type. However, it 
is possible to generalize our approach to arbitrary symmetries once the pressure 
dependence of the shear elastic constants (in addition to that of the longitudinal 
ones) is known. Further work in this direction is in progress. 
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